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ABSTRACT 

X-rays from planetary nebulae (PNs) are believed to originate from a shock driven 
into the fast stellar wind (v ~ 1000 km s _1 ) as it collides with an earlier circum- 
stellar slow wind (v ~ 10 km s" 1 ). In theory, the shocked fast wind (hot bubble) 
and the ambient cold nebula can remain separated by magnetic fields along a surface 
referred to as the contact discontinuity (CD) that inhibits diffusion and heat conduc- 
tion. The CD region is extremely difficult to probe directly owing to its small size 
and faint emission. This has largely left the study of CDs, stellar- shocks, and the 
associated micro-physics in the realm of theory. This paper presents spectroscopic 
evidence for ions from the hot bubble (kT « 100 eV) crossing the CD and penetrat- 
ing the cold nebular gas ( kT « 1 eV). Specifically, a narrow radiative recombination 
continuum (RRC) emission feature is identified in the high resolution X-ray spectrum 
of the PN BD+30°3639 indicating bare C VII ions are recombining with cool elec- 
trons at kT e = 1.7 ± 1.3 eV. An upper limit to the flux of the narrow RRC of H-like 
C VI is obtained as well. The RRCs are interpreted as due to C ions from the hot 
bubble of BD+30°3639 crossing the CD into the cold nebula, where they ultimately 
recombine with its cool electrons. The RRC flux ratio of C VII to C VI constrains the 
temperature jump across the CD to A kT > 80 eV, providing for the first time direct 
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evidence for the stark temperature disparity between the two sides of an astrophysical 
CD, and constraining the role of magnetic fields and heat conduction accordingly. Two 
colliding-wind binaries are noted to have similar RRCs suggesting a temperature jump 
and CD crossing by ions may be a common feature of stellar wind shocks. 

Subject headings: planetary nebulae: individual: BD+30°3639 - stars individual: BD+30 C 363 
- stars:winds - stars:Wolf-Rayet - X-rays:stars 


1. Introduction 

High velocity interactions between different astrophysical media due to collective plasma ef- 
fects on length scales much smaller than the mean free paths for particle collisions are commonly 
referred to as collisionless shocks. For extensive reviews see Tidman & Krall (1971); McKee & 
Hollenbach (1980). These shocks have been suggested to account for a wide variety of astrophys- 
ical phenomena ranging from the Earth’s weak bow shock to strong shocks in supernova remnants. 
Despite elaborate theoretical investigations, the detailed physics of collisionless shocks is far from 
being understood, mostly because direct observations or physical measurements of these shocks 
are notoriously difficult due to the extremely small size scales on which they occur. 

The extended X-ray emission in planetary nebulae (PNs) is most commonly explained by a 
(reverse) shock driven into the fast wind expelled by the central star as it collides with the earlier 
slow and massive wind (see review by Frank 1999). High angular resolution X-ray images of PNs 
seem to support this notion as the X-ray gas referred to generally as the hot bubble appears confined 
to the interior of the optical and IR nebula (e.g., Kastner et al. 2000). The low X-ray temperatures 
and observed luminosities cannot be explained by the present day wind velocities and mass loss 
rates. This suggests perhaps that the X-rays are due to the fast wind ejected during the early -PN 
phase or late post-asymptotic giant branch (post-AGB) phase, when the rapidly evolving wind was 
slower (v « 500 km s -1 ) and its mass loss rate was higher than it is today (Akashi et al. 2006). 
Alternatively, conduction of heat from the fast wind to the slow wind has been invoked to explain 
the observed low X-ray temperatures (Soker 1994; Steffen et al. 2008). In that case, it would be the 
nebular gas evaporating into the hot bubble that emits the X-rays. The fast PN wind running into 
the denser and slower nebular gas leads to relatively slow shocks (v f ~ 30 km s~') propagating in 
this dense gas. In a few cases, the morphology of the PNs suggest the fast wind is at least mildly 
collimated (Kastner et al. 2003; Sahai et al. 2003). 

High resolution X-ray spectroscopy of PNs is particularly challenging owing to the low X- 
ray flux from these sources. The only published X-ray grating observation of a PN is the Chan- 
r/ra/LETGS 300 ks exposure of BD+30°3639 (Yu et al. 2007, 2008). The analysis of Yu et al. 
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(2008) reveals a temperature range of (at least) 150-250 eV and extremely non-solar abundances. 
The elements C and Ne appear to be particularly enriched with respect to solar abundances (C/O 
30, Ne/O ~ 4), while Fe and N are deficient (Fe/O ~ 0.2, N/O ~ 0.4). In fact, these abundances 
found in the X-ray plasma closely trace the abundances measured directly from the wind of the 
present-day WC central star of BD+30°3639 indicating the X-rays do indeed originate from the 
fast wind, while evaporation of nebular-composition gas and possibly also conduction of heat to 
the surrounding nebula have no observable effect (Yu et al. 2008). 

The present paper is dedicated to an unusual emission feature in the LETGS spectrum of 
BD+30°3639 around 25.30 A, which we interpret as a narrow radiative recombination continuum 
(RRC) of bare C VII forming H-like C VI. This feature, observable only at the high spectral resolu- 
tion of the LETGS, provides a new insight into the region of the shocked fast wind of BD+30°3639 
and the micro-physics at play. 


2, Observations and Model 
2.1. The Spectrum 

BD+30°3639 was observed with Chandra for a total of ~ 300 ks, broken into several seg- 
ments from February - December 2006. All observations were carried out with the low energy 
transmission grating spectrometer (LETGS) and advanced CCD imaging spectrometer (ACIS-S) 
configuration. The data were processed using the standard Chandra pipeline and combined into 
one spectrum using the standard tools of the CIAO software package. The full details of the obser- 
vations and data reduction are given in Yu et al. (2008). 

As an initial step and following Yu et al. (2008), we employed the XSPEC package (version 
12.3.1, Amaud 1996) to fit the entire spectrum with two temperature components (2-T) of the astro- 
physical plasma emission code APEC (Smith et al. 2001), photoelectrically absorbed by cold gas 
using the absorption model of Morrison & McCammon (1983). A minor improvement over the Yu 
et al. (2008) approach is introduced for modeling the spectral line profiles; wavelength -dependent 
Gaussian smoothing is used. Since BD+30°3639 is ~ 4" across and spatially resolved by Chandra , 
and since LETGS is a dispersive slitless spectrometer, emission lines appear broadened in the spec- 
trum due to the target extent. This broadening is only slightly wider than the LETGS point-source 
line spread function and can be reasonably approximated by a Gaussian profile. In order to test 
the Gaussian spatial broadening approximation, we fit for the energy dependence of the smooth- 
ing width A E oc E a . The best-fit value obtained for a is 2.07 ± 0.15, consistent with spectral 
broadening of gratings in which AT is roughly independent of A, and thus A E oc AT/ A 2 oc E 2 . 

The best-fit spectral model is plotted over the source spectrum in Fig. 1. The most prominent 
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lines in the spectrum are those of C, O and Ne K-shell ions. Fe L-shell lines, often the best temper- 
ature indicators, are in fact not clearly identified in the spectrum due to the low Fe abundance. The 
best remaining temperature indicators are the line ratios of O VIII to O VII and of Ne X to Ne IX, 
which require here two non-degenerate temperature components of kT =160 eV and kT = 260 eV. 
The two temperatures likely represent a continuous temperature distribution in the hot bubble. For 
the case of a spherical wind, these range from the hottest inner region right behind the reverse 
shock, to the cooler, previously-heated, more extended regions near the contact discontinuity (CD) 
that are still much hotter than the few-eV outer nebula (see, e.g., Figs. 1 & 2 in Akashi et al. 2006). 
Much lower temperatures (kT <100 eV) are usually difficult to constrain from X-ray spectra, 
as the emissivities of most bright X-ray lines decreases strongly at these low temperatures. This 
problem is augmented here by the significant photoelectric absorption (Nh = 2.4 x 10 21 cm 2 ) to- 
wards BD+30°3639 and by the decreasing LETGS throughput with decreasing energy. The present 
spectrum, however, does provide unique temperature diagnostics as described in §2.3 and in §3.2. 

Relative abundances can be accurately measured with X-ray emission line spectra. The mea- 
surement of absolute abundances A Z /A H relative to hydrogen, on the other hand, is prone to sys- 
tematic uncertainties as it requires assumptions for the H contribution to the bremsstrahlung con- 
tinuum, which is harder to measure than narrow emission lines. In the case of BD+30 3639, the 
low temperature (low bremsstrahlung intensity), the low H abundance, and the generally non-solar 
composition of the X-ray plasma resulting from core He-buming, all make the assessment of A u 
even more uncertain. Consequently, we will limit our discussion in this paper to relative abun- 
dances, which are presented in Table 1. The present model parameters can be seen to agree with 
those of Yu et al. (2008) very well, as expected. For more details on the uncertainties associated 
with the absolute and relative abundance determination, see Yu et al. (2008). For the purpose of 
the present analysis, it is sufficient to note the unambiguously high C abundance in the hot bubble 
of BD+30°3639. 


2.2. The Carbon Footprint of BD+30°3639 

The 2-T model described above fits the overall spectrum quite well, but fails to reproduce 
the significant emission of what appears to be a dense complex of spectral lines in lire 25-26 A 
range. These lines coincide in wavelength with the high-order C VI Lyman series. The top panel 
in Fig. 2 zooms in on the relevant spectral region showing the data and the best-fit 2-T model, 
while highlighting the C VI Lyman series. The APEC model includes the C VI Lyman series, but 
only radiative decays from levels with n < 5, n being the principal quantum number of the upper 
level of the radiative transition. Thus, we first use the HULLAC code (Bar-Shalom et al. 2001 ) to 
complete the spectrum of the Lyman series up to n = 8 explicitly, and the decreasing contributions 
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of the remainder of the series by extrapolation. The fitted model with the complete Lvman series 
is plotted in the middle panel of Fig. 2. Evidently, the agreement of the model with the observed 
data is somewhat improved. 

Even after incorporating the complete Lyman series, the most prominent residual feature re- 
mains at 25.30 A. We identify this feature as the narrow radiative recombination continuum (RRC) 
of bare C VII recombining to form H-like C VI. In the following, we consider alternative identi- 
fications, but rule them all out. First, this feature need not be confused with the N VII Ly a line, 
which is very weak, but unambiguously resolved at 24.78 A (see Fig. 2). Indeed, Murashima et al. 
(2006), using a lower resolution spectrum of BD+30°3639, reported a high N/O abundance ratio. 
With the gratings, the two features are clearly distinct and the actual N/O ratio is much lower (Yu 
et al. 2008). Note that an overestimation of the N abundance should be expected from any low 
resolution observations that employ a model missing high-order C VI lines. Moreover, if a narrow 
RRC is present in the spectrum, it too would be confused at low resolution with N VII Ly a and 
would lead to a further overestimate of the N abundance. 

Other potential candidate lines around 25.30 A include the emission lines of L-shell ions of 
mid-Z elements such as Si, S, Ar, and Ca. The L-shell ions of S have no bright emission lines 
between 25 - 26 A (Lepson et al. 2003) and neither does Si (Lepson et al. in preparation). The 
strongest line of Li-like Ar XVI is at 25.02 A (Lepson et al. 2003), which given the LETGS 
resolving power, is easily distinguished from the observed emission feature at 25.30 A. Moreover, 
Ar XVI has equally bright lines at 24.87 and 23.53 A, neither of which is observed in the present 
spectrum. No other Ar ion has bright lines between 25 - 26 A. Our atomic computations show 
that Ne-like Ca XI has a 2p-4d line at 25.38 A. However, a Ca XI 2p-3d line at least five times 
as strong is predicted at 30.47 A, and should be accompanied by the bright 2p-3s lines at 35.65, 
35.67, and 35.79 A, none of which are observed in the spectrum. Given the extremely non-solar 
abundances observed in some PNs, we also checked the spectrum of Sc. O-like Sc XIV turns out 
to have a relatively bright 2p-3s line at 25.39 A. However again, stronger lines of Sc XIV that are 
not observed are expected based on our computations (e.g., a 2p-3d line at 22.95 A), as are lines 
of other L-shell Sc charge states. 

Charge exchange (CE) is another process that typically enhances high order lines when highly- 
charged ions mix with neutral atoms. In the context of astrophysical shocks, unshocked neutral 
atoms can easily penetrate the CD (electro)magnetic barrier and mix with the hot gas. If the 
enhancement is in sufficiently highly-excited levels of C VI (n > 10), the resulting lines could 
hypothetically fall below 25.5 A and would not be resolved from the C VII RRC in the LETGS 
spectrum. However, there are two strong arguments against such CE taking place in BD+30°3639. 
First, the UV radiation from the central star maintains high ionization (e.g., O III) in the inner parts 
of the cool nebular gas just outside the CD. Both BD+30°3639 and its sister PN NGC 40 feature 
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s trati fic ation in ionization, in which the innermost regions of the optical nebula, just outside the hot 
bubble, are the most highly ionized (Bryce & Mellema 1999; Sabbadin et al. 2000). Observations 
(and an estimate of the ionization parameter U « 0.03) imply that the fraction of neutral species 
available for CE with the X-ray ions is as small as < 1/1000. The second argument against CE in 
BD+30°3639 is that in order to produce an enhanced spectral signature below 25.5 A, the ioniza- 
tion potential of the neutral atom would need to be ridiculously low (< 2 eV, see Janev & Winter 
1985). This immediately rules out all of the abundant elements in BD+30°3639 including He, C, 
O, and Ne (as well as H) by a large margin and essentially all other elements as well. 

An intriguing possibility is CE with dust grains. Dust grains can survive UV charging and 
may have low ionization energies (i.e., work functions), which could be of the order of a few 
eV. Indeed, Matsumoto et al. (2008) find evidence for silicate grains in BD+30°3639 on scales of 
~ 4" that coincide with the diameter of the X-ray determined CD. However, at the corresponding 
distance of 4 x 10 16 cm from the central star, the grains inevitably are positively charged by more 
than a few elementary charges ( e ) as can be inferred from the work of Feuerbacher et al. (1973). 
This completely reduces the cross section for CE with positive (C VII) ions. In short, we deem the 
CE process as highly unlikely to produce the emission around 25.30 A. All of this leads us to the 
conclusion that the 25.30 A emission feature must be the RRC of C VII recombining to form C VI. 


2.3. The Carbon RRC 

The narrow width of the RRC, which is a direct measure of the electron energy distribution 
in the plasma, indicates recombination occurs with electrons of a few eV, two orders of magni- 
tude colder than those required to ionize carbon to its C VII state. In order to measure the exact 
recombining-electron temperature, we used the XSPEC redge model, which employs the flux den- 
sity profile 


Ie = F 


>RRC 


1 

kT c ° oli 


[photons s 1 cm 2 eV l ] 


( 1 ) 


where F RRC is the total photon flux in the RRC, kT c e °° l corresponds to the temperature in eV of 
recombining electrons, hv (= Ej+ electron energy) is the photon energy, and Ej is the ionization 
potential of the recombined ion. The ionization energy of C VI was set at £/ = 490.02 eV and the 
two remaining independent parameters were fitted, yielding kT c e °° l - 1.7 ± 1.3 eV and a total RRC 
flux of F RRC n = (8.8 ± 3)xl0“ 5 photons s" 1 cm -2 . The upper limit on kT c e °° l is robust as the maximum 
width of the RRC is tightly constrained by the data. Conversely, the lower limit is not as well 
determined due to the RRC being just broader than the instrumental line spread function expected 
from the extended angular size of BD+30°3639. At the distance of BD+30°3639, d - 1 2 kpc (Li et 
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al. 2002), we find the isotropic rate of C VII recombination events resulting in photons of - 25.30 A 
(i.e., with cool electrons) to be 

lev,, = 4 nd 2 F R c R v c „ = (1.5 + 0.5) x 10 40 s -1 . ( 2 ) 

For physical consistency, we also include in the model the C VI line intensities due to recombi- 
nation onto excited levels (up to n = 8) and ensuing radiative cascades, with relative contributions 
that are appropriate for kT c °° l . In the recombination model for kT c e °° l =1.7 eV, the unabsorbed 
relative intensities of the RRC, Lya, Ly/?, Lyy, and Ly 8 are 1.0, 1.35, 0.22, 0.08, and 0.04, re- 
spectively. These contributions were calculated using HULLAC (Bar-Shalom et al. 2001 ). and 
here added to the model as narrow Gaussians at their fixed wavelengths. The final model around 
the relevant spectral region can be seen in the lower panel of Fig. 2, including the RRC and the 
low contribution of recombination to the lines compared to the high-7 emission. Note that photo- 
electric absorption in this part of the spectrum is significant (N H = 2.4 X 10 21 cm -2 ), and decreases 
strongly away from the C I edge at ~ 43.5 A towards shorter wavelengths, so that low order lines 
of the series (Ly a, Ly/?) are much more absorbed than the high order lines. The best-fit flux model 
corrected for interstellar absorption is presented in Fig. 3 and the model parameters are listed in 
Table 1. 

The miniscule flux due to recombination of highly ionized atoms with cool electrons com- 
pared with the strong line emission originating from hot gas makes the detection of the former 
tremendously hard. It is obvious from Figs. 2 and 3 that if not for the unusually high C abundance 
( A c /A 0 « 40) in BD+30°3639, this detection would not have been possible. Nevertheless, we did 
look for more evidence for recombination in the spectrum. The RRC of C VI forming C V falls at 
31.63 A. Due to the low count rate at these long wavelengths, only an upper limit can be obtained, 
which is: 3.7 X 10 -4 photons s -1 cm -2 (90% confidence). This flux is high (not very constraining), 
particularly when compared with the C VII RRC flux of (8.8 ± 3) x 10 -5 photons s ' cm 2 , The 
difficulty is due to the low effective area of LETGS at this wavelength (~ 2 cm 2 ) and the strong ab- 
sorption near the C I edge. Nonetheless, this upper limit still proves useful, as will be demonstrated 
in §3. The flux in the O VII forbidden line at 22.01 A and the He/? line at 18.63 A can be seen 
in Fig. 1 to exceed the flux predicted by the hot plasma model, which constitutes further, though 
tentative, evidence for recombination as the upper levels of these lines are preferentially populated 
by recombination and radiative cascades. The O VIII RRC forming O VII falls at 16.78 A. but 
only an upper limit (90% confidence) to the flux of 7.6 X 10 -6 photons s -1 cm -2 can be obtained. 
Blending with the Fe XVII line at that wavelength is possible, but not expected to be significant, 
since Fe is considerably under-abundant. Although unambiguous detection is possible here only 
for the highly overabundant C, the upper limits are all consistent with the coolest (kT - 100 eV) 
species in the hot bubble recombining with kT »2eV electrons in the nebular gas. 
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The population of intermediate-temperature electrons ( kT > 3 eV) cannot be significant, as it 
would have smeared out the high contrast of the C VII RRC. To demonstrate this, we attempted to 
fit the same model as before, but with two C VII RRC components instead of one. We note that the 
quality of the data does not warrant a second RRC (see Fig. 2). Nevertheless, the main component 
was fixed at its previous best- fit value of kT\ = 1.7 eV, while another cool component kT 2 was 
tested for (not to be confused with the high temperatures of the 2-T plasma model of §2.1). The 
flux of both components was allowed to vary in the fit, as were the other model parameters. In 
Fig. 4, we plot lcr flux limits of each of the two RRCs as a function of kT 2 . The 1.7 eV component 
is significant at all kT 2 values, corroborating the existence of cool electrons. At high temperatures 
0 kT 2 > 20 eV), the solution gradually converges to the single RRC result. A 10 - 20 eV RRC 
is possible, although it is consistently weaker than the 1.7 eV RRC (Fig. 4). Below ~i0 eV, the 
plot shows that there is a high degree of confusion between the two components. If we allow both 
RRC temperatures to vary, a best fit is obtained for kT , = 1 eV, kT 2 = 10 eV, I\ = 1^, / 2 = I 
10 40 recombination-photons s -1 . These do not represent distinct components as the temperatures 
are highly confused and interchangeable. Only an upper limit of kT < 20 eV can be obtained 
for the hotter components. The addition of the second RRC (two extra free parameters) lowers 
the reduced-^ 2 from 0.782 to 0.775 (A^ 2 = 1.97), if the method of Gehrels (1986) is used. We 
conclude that recombination with electrons hotter than ~ 3 eV is not required by the data (see 
Fig. 2). However, insignificant recombination with electrons of up to ~ 20 eV cannot be ruled out. 
A good probe of these intermediate temperatures could be O VI emission (Gruendl et al. 2004), 
but we are not aware of O VI measurements for BD+30°3639. 


3. Discussion 

In order to produce the observed narrow RRCs, highly-charged ions have to collide with cool 
electrons at temperatures corresponding to only a few eV. Narrow RRCs are typical of photo- 
ionized X-ray plasmas, where the ionization level is high, but the electron temperature remains 
low; good examples are X-ray binaries, or Seyfert 2 ionization cones. In contrast to such sources, 
PNs have no hard enough X-ray source to ionize carbon up to the C VII state. One might think 
that plasma instability around shocks could entail rapid radiative cooling in which ionized atoms 
recombine with cool electrons. However, recombination is one of the primary cooling mechanisms, 
and highly charged species such as C VII would have long recombined by the time the plasma 
cools from above 100 eV to a few eV (Sutherland & Dopita 1993, Fig. 16 therein), which rules 
out instability and rapid cooling as the origin of the RRCs. In principle, adiabatic expansion could 
produce RRCs if it is sufficiently rapid that ions do not have time to recombine before the plasma 
cools to a few eV, i.e., ’’ionization freezing”. In reality, this is unlikely. The expansion time depends 
on the geometry, but the confined geometry of the hot bubble by the dense nebula does not favor 
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such rapid expansion. During the expansion, due to the adiabatic relation T oc /? 2 \ and since 
recombination rates (Eq. 6 below) scale as nT~ l/2 oc T, recombination would preferentially occur 
during the hot, dense phases and not during the low-density phases when the plasma cools down to 
a few eV. Hence, it would not produce the observed narrow RRCs. Moreover, for the hot gas to cool 
from 100 eV to 1 eV, its density would have to decrease by a factor of 1000. Consequently, in order 
to produce the observed emission measure of the RRC (n 2 V ~ 10 53 cm -3 ), the size of the source 
would have to be more than two orders of magnitude larger than the hot bubble, which can be 
ruled out. A transient ionizing scenario, as expected for shocked gas, cannot produce RRCs either, 
although the initial ion temperature would be much higher than the electron temperature. This is 
because electron impact ionization follows electron heating. Even if one invokes a yet-unspecified 
mechanism to ionize the hot atoms (up to C VII) before the electrons heat up, the collisions of cold 
electrons with hot ions much more efficiently heat the electrons than result in recombination. We 
conclude that the most plausible physical origin of the observed RRCs is highly-charged species 
that are heated and ionized in the hot bubble (kT >100 eV) and then interact with much cooler 
electrons (kT ~ a few eV) in a distinct medium. 


3.1. Ions Crossing The Contact Discontinuity 

The most obvious reservoir of cool electrons lies in the cold nebula, where typical temper- 
atures indeed correspond to a few eV. The nebula is photo-ionized by the central UV source, so 
unbound electrons are abundant. The observed recombination can take place, therefore, if ions 
from the hot bubble cross the CD plane and penetrate the slow wind where the cool electrons re- 
side. In the simplified hot-bubble picture, the outer part of the hot bubble is also its densest and 
relatively coolest region (for the spherically symmetric case see, e.g., Akashi et al. 2006 , 2007 ). 
This picture is consistent with the lowest charge states in the hot bubble crossing the CD and re- 
combining, as the spectrum seems to suggest. The narrow RRC and the consequential limits on 
intermediate temperatures should help constrain the role of heat conduction by electrons as well as 
the role of nebular evaporation into the hot bubble. 

Heat conduction by electrons can be suppressed even by magnetic fields as weak as 0. lpG 
(Soker 1994). The viability of such magnetic fields has not been established, but if they do exist, 
they could preserve the steep temperature gradient across the CD required to produce a narrow 
RRC. The length scales for such a gradient are of the order of the electron Larmor radius, while the 
Larmor radius for ions can be much larger. For particles with mass m and charge q at a temperature 
T, Rl ~ V !q. Thus, ions can effectively penetrate the cold nebula across a magnetically held 
CD, while the opposite effect of cold electrons penetrating the hot bubble would be negligible. 
The Larmor radius of a bare C VII ion with kT hot « 100 eV gyrating along a typical interstellar 
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magnetic field component of B = 1 jjG (parallel to the CD surface) is 


K l=6X io»f^n— I 1 ' 2 ™ 

\lj uG) \l00eV / 


The corresponding half-circle (cross-back) rotation time would be: 


ncM c ( B N 1 

t l = — — ■ ■ = 650 1 - — - | s 


6 eB 


1/iG 


( 3 ) 


( 4 ) 


where M c is the atomic mass of carbon. For comparison, the slow-down time by Coulomb colli- 
sions of a C VII ion with the ambient cold gas can be estimated from Eq. (5-29) of Spitzer ( 1 956) 
to be (independent of ion temperature) roughly: 


25 


cool 


k T co°l\ V2 I 

leV j \10 4 cm 


/ m cool\ l/ 2 




(5) 


where m co ° l and n co ° l are the mass and number density, respectively, of the cool-gas species respon- 
sible for slowing down the fast C ion. More precisely, for the temperature and density of Eq. (5), if 
the C ion collides with cool C atoms r* « 4 s, with He atoms t s « 10 s, and the longest it may take 
to slow down, if it is stopped by electrons, is r s « 1200 s. Obviously, the strength of the magnetic 
field and to a lesser extent the density in Eqs. (4, 5) are fairly uncertain (even to an order of magni- 
tude), but the fact that r s is less than or comparable to r L implies that ions crossing the CD have a 
good chance of being stopped by the cool plasma, rather than returning to the hot bubble. Chances 
can be much higher, if the magnetic fields are weaker than 1 jjG, or absent altogether. Note that in 
order for the hot C ions to reach the 1 .7 eV gas, they need to cross any intermediate region (say 
kT = 10 eV) without being stopped. This considerably limits the size of such a region through 
t s (T ) and the ion velocity, which for intermediate temperature kT = 10 eV yields ~ 1 0 9 cm. 


Compared to these short time scales, the recombination of C VII with the cool nebular elec- 
trons can take several months: 


T rec 


1 


1.8 x 10 7 


kT: 


cool \ 1/2 


leV 


) (ic^r 


(6) 


where = 5.5 x 1 0~ l2 cm 3 s~ 1 (k7 e c " o 7l eV)~ 1/2 is the radiative recombination rate coefficient. 
These time scales suggest the following plausible scenario: Ions from the hot bubble (k T 1 -100 e V) 
cross the CD into the cool nebula. Many of them slow down before they can gyrate back, and ther- 
malize by collisions with the nebular electrons, cooling down to a few eV. After a much longer 
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time, and perhaps after they diffuse further upstream, the ions eventually recombine with the cooi 
nebular electrons. 


3.2. Temperature Diagnostics at the Contact Discontiuity 

As seen above, the width of the RRC can provide a tight constraint on the electron temperature 
on the nebular side of the CD. It would also be useful to measure the temperature of the hot plasma 
close to the CD, but just on its hot side. This temperature can be readily obtained from the limit on 
the recombination rate of C VI. In a steady state, the rate of recombination with cool electrons is 
balanced by the ion crossing rate. Indeed, steady state is reached within the typical recombination 
times T rec < year (Eq. 6), which are much shorter than the PN age. The measured recombination 
rates I rec , therefore, reflect the corresponding ion densities n l0n and thus the fractional abundances 
in the hot gas fi 0n (T hot ) °c n wn . The C VII to C VI ionic density ratio is particularly sensitive to 
the temperature around kT ~ 100 eV. These densities and temperatures refer to the edge of the 
hot bubble, right at the CD surface. Taking into account that C VII ions crossing the CD will 
produce two photons, one in the C VII RRC and subsequently one in the C VI RRC, one can write 
Ir cv„ x fcv„(T hot ) and 7™, oc fcv,(T h "‘) + fcv„(.T ho '), or: 

/ / Thot\ Tree 

JCVIlU ) _ l CVII 

fcv,(T hot ) ~ I r c %, - I r c ec vn 

The measured upper limit for Vg c vl (§2.2) places a lower limit on the ratio fcvttij hat )l fc vi(T h " 1 ) > 
0.31 ± 0.1. The theoretical ionic fraction ratio as a function of temperature (Mazzotta et al. 1998) 
is plotted in Fig. 5 along with the measured lower limit imposed by Eq. (7). These results indicate 
that the temperature on the immediate hot side of the CD must be at least kT ho ‘ > 88 + 4 eV. 
Together with the accurate measurement of kT C0 ° l = 1.7 ± 1.3 eV from the width of the RRC, 
this implies a temperature jump corresponding to at least 80 eV across the discontinuity, with no 
significant plasma at intermediate temperatures. In principle, the O RRCs could have provided 
independent diagnostics, perhaps even a valuable upper limit on T hot , since the A c /A 0 abundance 
ratio in the hot bubble is well constrained. Unfortunately, the present upper limit on the O VIII 
recombination is not restrictive enough to provide meaningful constraints. 

The requirement for a steep temperature jump corresponding to ~ 80 eV across the CD is 
robust and it raises an intriguing question regarding the role of magnetic fields and heat conduction. 
Heat conduction would drastically change the temperature and density profiles across the CD. 
Instead of a sharp-contrast interface, a more gradual conduction front would manifest the transition 
between the hot bubble and the nebula as nebular gas evaporates into the hot bubble. On the face of 
it, this picture seems to be in contrast with the LETGS spectrum of BD+30°3639 and the present 
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findings. In addition to the sharp temperature jump across the CD implied by the present analysis, 
the fast-wind composition of the hot bubble in BD+30°3639 argues against a significant role for 
nebular evaporation into the hot bubble, at least in this source (Yu et al. 2008). Note, however, that 
the penetration of hot-bubble material into the nebula will alter the chemical abundances of the 
immediate inner parts of the nebula, eventually changing the composition to that of the fast stellar 
wind. Such abundance gradients in the inner nebular regions, gradually changing from fast-wind 
composition to nebular composition as the distance from the center increases, were indeed found 
by Sabbadin et al. (2000) in spatially resolved observations of NGC 40. 

Steffen et al. (2008) published numerical simulations for PN hot bubbles that include heat con- 
duction in the magnetic-free limit. Strictly speaking, these models cannot be applied to BD-f-30 3639. 
since they assume solar abundances dominated by hydrogen, while the X-ray gas in BD+30 3639 
considerably departs from solar composition, consists predominantly of He and C, and is extremely 
H-deficient. Nevertheless, we note that the temperature profiles of Steffen et al. (2008) at the outer 
edge of the conduction front do appear to sustain a sharp temperature drop by approximately two 
orders of magnitude (see their Fig. 4), perhaps due to the rapid cooling of plasma between a few 
eV and 100 eV. However, the spatial resolution of these simulations is of the order of 10 15 cm, 
while the thickness of the transition layer estimated in the present work is several times 10 13 cm 
(§3.3). Careful modeling of heat conduction on the small scales probed here and under conditions 
closer to those of BD+30°3639 (He and C dominated plasma) will be needed before making a more 
conclusive and quantitative statement regarding the role of magnetic fields and heat conduction in 
PNs with WC central stars. 


3.3. Further Estimates: Carbon Mass Outflow, Densities, and Penetration Depth 

The measured RRC intensity (Eq. 2) implies a CD mass crossing rate of M C vn = (4.7 ± 1 .6) x 
10 -9 M 0 yr -1 , and further dividing by the fractional abundance of C VII at 100 eV (fcvn= 0.27, 
Mazzotta et al. 1998), yields a total recombining C mass rate of Me * (1.8 ± 0.6) X 10 '' M 0 yr 1 . 
We recall that approximately 50% of the hot bubble mass of BD+30°3639 is likely in the form C 
(Yu et al. 2008; Marcolino et al. 2007). If this estimate holds near the CD, it implies a mass crossing 
rate of M « (3.5 ± 1.2) X 10 -8 M 0 yr -1 . The most recent fast-wind mass outflow rate estimate for 
BD+30°3639 is 5 x 10 -7 M 0 yr -1 (Marcolino et al. 2007), which is an order of magnitude less than 
the estimate of Leuenhagen et al. (1996). Therefore, the CD mass crossing rate is, at the most, less 
than 10% of the mass outflow rate and even a smaller fraction of the early PN wind that had an 
even higher M. To that end, the ions crossing the CD seem not to have an appreciable effect on the 
stellar wind, or the shock dynamics. 

With a simplifying assumption of a smooth CD geometry that ignores wind clumping (c.f., 
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Prinja et al. 2007), the ion density can be inferred from the observed recombination rate. We 
assume that ions cross from the hot bubble across the CD into the cold nebula exclusively due to 
their thermal motion. As in the previous section, in a steady state in which the C VII crossing rate 
is balanced by the observed recombination rate one can write: 


I<cvn ~ 0.5v ± S cd^cvh ( 8 ) 

where v ± = ^kT hot /M c is the average magnitude of the C VII velocity component perpendicular 
to the CD, whose surface area is Scd • The factor of 0.5 comes from the isotropy of the thermal 
motion. Eq. (8) can then be used to express the C VII ionic density on the hot side of the CD as 
follows: 


n 


hot 

CVII 


x5 ( 


'CD 


2x 10 34 cm 2 ) 


kT hot X-l/2 

100 eV/ Cm 


(9) 


where we assumed the CD surface area Scd to be a perfect sphere with the observed X-ray radius 
of BD+30°3639: Rcd = 4 X 10 16 cm (2" at 1.2 kpc). Using again the ionic fraction fcvi/i 100 
eV) = 0.27 yields a C density of nc ~ 2 cm -3 . Assuming as before and following Marcolino et al. 
(2007) and Yu et al. (2008), a non-solar composition of gas depleted of hydrogen and dominated 
by He and C (Ac/Ane ~ 0.4 by number), we can estimate the electron density in the hot bubble 
to be n e « 3 n He « 8n c ~ 16 cm -3 . To that extent, the observed recombination rates appear to- 
be consistent with the scenario of shocked hot bubble gas penetrating the slow wind, and with the 
theoretical models and numerical simulations of such a scenario. 


Of course, the assumption of spherical symmetry in these models and in Eq. (9) is proba- 
bly an oversimplification. Collimated fast winds would imply higher densities. Furthermore, the 
simplified picture that only includes crossing of hot ions to the cool side must be incomplete as 
it would breach the plasma neutrality and would very quickly induce a strong electric field that 
would inhibit further ion crossings. Hence, there must be a reverse-charge current, for example 
of electrons following the ions, or nebular ions entering the hot bubble. Further discussion of the 
complex morphology of the electric and magnetic dynamical structure around the CD is clearly 
beyond the scope of our study, but these findings could be useful to constrain such investigations 
in the future. 

Finally, the C VII RRC intensity, along with the above estimate of the density, can be used 
to approximate the penetration depth A R of the hot-bubble material into the nebula. Using the 
fact that the recombination rate /£“„ scales with the electron and ion densities integrated over the 
entire emitting volume, which in turn can be approximated by a thick sphere (S C d) of width A R. 
just outside the observed X-ray hot bubble (the CD), one can write 
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A R = 


Tree 

_CV77 

n e a RR S cd (ncvn) co ° l 


( 10 ) 


where ( n C vii) co ° 1 is the steady-state average C VII density on the cold side of the CD. This average 
likely represents a density gradient that cannot be constrained without a valid model for diffusion 
and plasma effects in the nebula. All we can do at this point is scale A R with (n C vii) co "', which can 
be arbitrarily equated with n h ^‘ vu of Eq. 9, to yield 



At the distance of BD+30°3639, this is only ~ lmas on the sky, which clearly cannot be resolved 
by X-ray telescopes even if the actual value of (n C vii) co ° l is much smaller. Such a shell needs to 
cool rapidly enough to radiate away the heat due to the penetration of the hot gas. A penetration 
rate of ~ 10 42 particles per second (c.f., Eq. 2), each contributing 100 eV, yields a heating rate of 
~ 10 32 erg s' 1 . If the density in the nebula is ~ 10 4 cm -3 , a shell at a temperature of kT ~ 2 eV 
and with a volume of S C d&R ~ 6x 10 47 cm 3 would cool at a rate of at least 6xl0 33 erg s' 1 (given 
A > 10' 22 erg cm 3 s'\ Sutherland & Dopita 1993), safely above the heating rate. Turning the 
argument around and requiring that the shell be able expel the heat from penetrating ions constrains 
its thickness to a more conservative value of AR > 5 x 10 11 cm. 

There is reason to believe the stellar wind material continues to diffuse into the nebula much 
after it recombines as spatially resolved images of NGC 40 show abundance gradients in low- 
ionization species on scales of a few arcseconds, or a few times 10 16 cm. The penetration of the 
stellar wind from the hot bubble deep into the nebula is a different effect from nebular evapora- 
tion, although both processes similarly defy the simplified notion of stark disparity in composition 
between the two media. In contrast with the chemical composition, however, the temperature 
disparity does appear to preserve the discontinuous nature of the interface. 


3.4. Ion Crossing of the Contact Discontinuity as a 
General Phenomenon in Stellar- Wind Shocks 

The present discovery of a narrow carbon RRC in a PN shock is possible owing to the cold, 
yet ionized gas of the ambient nebula, but also to the extremely high C abundance in the stellar 
wind and the hot bubble. It would be natural to suppose, however, that hot ions crossing the CD 
should be a more general phenomena of stellar-wind shocks. Indeed, we found two more cases 
of narrow RRCs in hot plasma sources; both are colliding wind binaries with a C-rich Wolf-Rayet 
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star and an O star companion: y 2 Velorum and 9 Muscae. 

y 2 Velorum is a WC8+07.5 stellar binary system, which like BD+30°3639 (WC9) is highly 
enriched in C, O, and Ne. Schild et al. (2004) report the detection of C VII and C VI RRCs in the X- 
ray spectrum of y 2 Velorum, and upper limits to the RRCs of oxygen ions. Similar to BD+30°3639, 
the observed recombining gas in y 2 Velorum reflects the WR chemical composition. In the binary 
case, the O star wind plays the role the cold nebula plays in PNs, essentially stopping the WR wind 
and sending back a reverse shock. The electron temperature of the C VII RRC measured by Schild 
et al. (2004) is 3.3 + 0.7 eV, consistent with that of an un-shocked O star or WR wind. The RRCs 
in y 2 Velorum persist when the thick WR wind absorbs the inner, hot X-ray components implying 
they originate downstream from the stagnation point along the shock front (e.g., region 3 in Fig. 9 
of Schild et al. 2004). Schild et al. (2004) speculate that the narrow RRCs could be due to cold 
plasma in extended regions of the wind, photo-ionized by hard X-rays from the hot wind-collision 
region. However, they also note significant problems in finding a model that can account for both 
the high ionization parameter required to photo-ionize C to its highest charge state (suggestive of 
low density) and the high emission measure of the recombining plasma (high density). 

Schild et al. (2004) also raise the possibility that recombination may occur following adiabatic 
cooling of the shocked gas. As discussed in §3, this would require the (adiabatic) dynamical times 
to be considerably shorter than the recombination times, which is highly unlikely in a dense stellar 
wind. Using the parameters of Schild et al. (2004), we estimate the density 1 AU away from the 
WR star to be n e « 4 X 10 8 cm" 3 , which implies C VII recombination times of r m « 5000 s (for 
kT = 100 eV). Adiabatic expansion times are harder to estimate without knowing the geometry, 
but the dynamical time it takes plasma that cools down to kT =1 eV ( C soun d ~ 10 km s” 1 ) to 
expand to say ~1 AU is roughly 10 7 s » T rec . As in the PN case discussed in §3, if the hot gas 
expands and cools from kT > 100 eV to a few eV, the density would have decreased by a factor 
of ~ 1000. This would require an enormous volume to produce the observed RRCs. We suggest 
alternatively that the RRCs in y 2 Velorum, in analogy with the case of BD+30°3639, are due to 
highly-ionized C atoms that were shocked in the WR wind, crossed the CD surface, and interacted 
with un-shocked electrons of the much cooler, yet also ionized plasma from the O star wind. The 
high carbon abundance of the WR wind, as opposed to the much lower abundance in the O star 
wind, suggests that the hot ions must come from the WR side. The wind from the O star is also 
shocked (and hot) close to the stagnation point and therefore this hot-WR-ions/cold-O-electrons 
interface must exist farther away along the bow as Schild et al. (2004) infer from observations. 

Sugawara et al. (2008) identified narrow C VII and O VIII RRCs in the X-ray spectrum of 
another colliding wind binary, 6 Muscae, a WC6+09.5 system fairly similar to y 2 Velorum. They 
suggest the possibility of ions escaping the bow shock layer and interacting with the wind, but 
provide no further discussion. Again, we suggest the cool recombining gas interacting with the 
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highly-ionized plasma in 6 Muscae be interpreted as evidence for shock heated plasma in the metal- 
rich WR wind crossing the CD and mixing with the cooler O star wind. The bright RRCs in the two 
stellar binary systems and in BD+30°3639 suggest this kind of mixing and perhaps other micro- 
physical processes around the CD are common to stellar-wind shocks in different systems ranging 
in size from ~ 10 13 cm in stellar binaries to 10 16 - 10 17 cm in PNs. The high C abundance in these 
sources obviously facilitates the detection of the RRCs of C VII and C VI, which are relatively 
isolated in the spectrum. O RRCs are harder to identify unambiguously owing to blending with 
Fe L-shell lines. We suspect, nonetheless, that the observed RRCs are traces of shocked plasma 
crossing the CD irrespective of element biases. 

An interesting question is whether similar CD crossing takes place in young sheil-like super- 
nova remnants (SNRs) and whether it leaves a measurable signature that can give an idea of the 
temperature jump across the CD, as it does in BD+30°3639 and in the colliding binary winds. The 
fast stellar wind colliding with the PN is qualitatively analogous to the SN ejecta slamming into 
the circum-stellar material in young SNRs. However, due to the high velocity of the SN explosion 
and the low density of the circum-stellar gas, temperatures on both sides of the CD, heated by 
the forward and reverse shocks, are both expected to be high (kT > 100 eV and possibly much 
higher). Additionally, the immediate region around the CD represents the oldest shocked plasma, 
so the electrons that were initially cool have had time to equilibrate with the hotter protons. Con- 
sequently, ions crossing the CD would not produce a narrow RRC even if they recombine. More- 
over, unlike in PNs, the influence of strong magnetic fields on the CD in SNRs is not in doubt, 
as indicated by radio synchrotron emission from thin filaments attributed to the CD region. Tire 
interstellar magnetic field, and in particular its tangential component, can be enhanced by a factor 
of a few at the CD due to field compression by simple charge advection (Cassam-Chenai et al. 
2005), or by up to two orders of magnitude due to hydrodynamical instabilities (Jun et al. 1995). 
These fields are more than sufficient to sustain the temperature jump and perhaps even to totally 
suppress ion crossing. By substituting a magnetic field of 100 pG and a density of 1 cm" 3 into 
Eqs. (4, 5), one easily sees that a thermal ion gyrating across the CD of a SNR will spend far 
too little time on the cool side of the CD (t l ~ seconds) in order for it to significantly cool down 
and remain on the other side (r s ~ days), even if cool electrons were present. Indeed, the highest 
quality published X-ray grating spectra (especially beyond 25 A) of young SNRs do not appear to 
feature detectable narrow RRC emission (Rasmussen et al. 2001; Behar et al. 2001). We therefore 
conclude that narrow RRCs due to ion crossing of the CD are likely a unique signature of relatively 
slow shocks in dense stellar wind sources. 
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4. Conclusions 

BD+30°3639 was observed with Chandra LETGS for a total of 300 ks providing the first and 
so far only high-resolution X-ray spectrum of a PN. As already shown by Yu et al. (2008), the 
spectrum can be fitted by an absorbed, two-temperature plasma model, though it likely represents 
a distribution of temperatures between ~100 - 300 eV. The spectrum is dominated by C, Ne, and O 
emission lines with very little N and Fe. The Carbon abundance is particularly high, indicating the 
X-ray gas originates primarily from the present-day WC stellar-wind and making the C emission 
features most suitable for detailed analysis. 

We detect in the spectrum a narrow RRC of bare C VII forming C VI by recombination with 
cool ( kT e = 1.7 ± 1.3 eV) electrons. We suspect a bright RRC of H-like C VI forming C V also 
appears in the spectrum, but the low S/N of that feature allows only an upper limit to its flux. We 
interpret the RRCs as direct evidence for penetration of hot, highly ionized plasma into the cool 
nebula past the contact discontinuity, which separates the heated fast wind from the slow wind. 
This finding requires a steep temperature gradient, as intermediate electron temperatures higher 
than kT e — 3 eV would have broadened the RRC beyond its observed width. Using the measured 
lower limit to the flux ratio of the C VII and C VI RRCs, the temperature difference between 
the hot and cold plasma across the interface is found to be A kT > 80 eV. Such a steep gradient 
can be sustained by magnetic fields, which would significantly suppress heat conduction. On the 
other hand, magnetic fields may not be needed to preserve the sharp temperaure drop between 
the hot bubble and the nebula. Numerical simulations of heat conduction in the literature (Steffen 
et al. 2008) do appear to preserve a sharp gradient, though their grid resolution is too low for the 
processes discussed above. Simulations with the appropriate WC wind composition and at higher 
resolutions will need to be confronted with the present results in order to quantitatively test the role 
of magnetic fields in PNs. 

The measured recombination rate of (1.5 ± 0.5)xl0 40 s _I implies a mass crossing rate of 
~ 3.5 xl0~ 8 M 0 yr -1 , which can be a few percent of the total fast-wind mass, implying the fast wind 
gas is successfully making its way through the nebula despite the putative magnetic barrier. The 
deep penetration of stellar wind gas into the optical nebula is indirectly supported by abundance 
gradients observed in the innermost regions of NGC 40 (Sabbadin et al. 2000). This picture is quite 
different from that of nebular evaporation into the hot bubble for which we find no observable sign. 
Finally, we note two stellar (WC + O) binary X-ray sources with similar RRC features. We take 
this as evidence for the generality of ions crossing CD surfaces in slow, dense stellar wind shocks, 
regardless of whether the shock front is on the scales of a massive-star binary system or a PN. 
The hotter and more highly magnetized environments of CDs in supernova remnants, on the other 
hand, are less favorable for ion crossing, or for producing narrow RRCs. 
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Fig. 1.- — LETGS background-subtracted spectrum of BD+30°3639, using bin size of 0.1 A, over- 
laid with the best-fitted model 
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Table 1 : Parameters of the best-fitted model 


Parameter 

N// 

kT\ 

EM\ 

kT 2 

EM 2 

Ac/Ao 

An/Ao 

Ao 

ANe/Ao 
Ams/Ao 
Asi/Ao 
Af e j Ao 


Value ± Error Yu et al. (2008) 11 
0.24 ±0,04 | 0,24 ±0.04 

Thermal Components 
260 ± 30 ^ 250 ± 30 I 


260 ± 30 
50 ± 10 
160 ± 30 
160 ± 45 
37 ±8 
0.3 ± 0.4 
0.9 ± 0.4 
3.8 ±0.4 
0.8 ±0.5 
2.0 ± 1.1 
0.15 ±0.1 


70 ±50 
150 ±30 
170 ± 35 
33 ± 15 
0.4 ± 0.4 
0.9 ± 0.4 
3.8 ± 0.8 
0.7 ± 0.5 

0.2 ± 0.2 


Units 

10 22 cm -2 

eV 

10 53 cm -3 
eV 

10 53 cm -3 
solar 1 ' 
solar 1 ' 
solar 1 
solar 1 
solar 1 
solar 1 
solar 1 


Recombination Component 


t Anders & Grevesse (1989) 

NOTE: Abundance uncertainties are 90% confidence range 




kjcool 

1.7 ±1.3 

eV 

f'RRC 
r CVII 

8.8 ±3 

10 -5 photons s _1 cm -2 
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Fig. 2. — A zoom into the Carbon spectral region with 0.1 A bin resolution. Top: Best-fit 2-T 
APEC model. Middle: 2-T APEC model supplemented with the complete C VI Lyman series. 
Note the inadequacy of this model to fit the feature around the C VII RRC at 25.30 A. Bottom: The 
full model (red) including emission following recombination with cool electrons (thinner brown 
line). Note that interstellar absorption increases drastically with wavelength across this range, 
which explains the skewed Lyman series line ratios. See Fig. 3 for the fluxed model. 
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Fig. 3— Fluxed best-fit spectral model forBD+30°3639. The small contributions of C VII recom- 
bination with cool electrons (as well as the upper limit for the C VI RRC - dashed line) compared 
with the strong line emission originating predominantly from hot gas demonstrate how hard it 
is to detect recombination from ions crossing CD surfaces in astrophysical shocks and how the 
extremely high C abundance of BD+30°3639 (A c /A 0 « 40) is crucial for the present detection. 
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Fig- 4. — Constraints on C VII recombination rates from two electon-temperature components as a 
function of kT 2 the temperature of the second component. The temperature of the first component 
is held constant at kT = 1 .7 eV. 







